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Abstract. The nature of the dark components (dark matter and dark energy) that dominate the current cosmic evolution is a 
completely open question at present. In reality, we do not even know if they really constitute two separated substances. In this 
. paper we use the recent Cosmic All Sky Survey (CLASS) tensing sample to test the predictions of one of the candidates for 

' a unified dark matter/energy scenario, the so-called generalized Chaplygin gas (Cg) which is parametrized by an equation of 

I state p = -A/p'^^ where A and a are arbitrary constants. We show that, although the model is in good agreement with this radio 

source gravitational lensing sample, the limits obtained from CLASS statistics are only marginally compatible with the ones 
I ' obtained from other cosmological tests. We also investigate the constraints on the free parameters of the model from a joint 

O I analysis between CLASS and supernova data. 

'• 

^ ' Key words. Cosmology: theory - dark matter - cosmological parameters 

1. Introduction 

^ As is well known, there is mounting observational evidence that our Universe is presently dominated by two exotic forms of 
■ " " matter or energy. Cold, nonbaryonic dark matter, which accounts for ^ 30% of the critical mass density and whose leading 
particle candidates are the axions and the neutralinos, was originally proposed to explain the general behavior of galactic rotation 
curves that differ significantly from the one predicted by Newtonian mechanics. Later on, it was also realized that the same 
concept is necessary for explaining the evolution of the observed structure in the Universe from density inhomogeneities of 
the size detected by a number of Cosmic Microwave Background (CMB) experiments. Dark energy or quintessence, which 
accounts for ^ 70% of the critical mass density and whose leading candidates are a cosmological constant A and a relic scalar 
field (p, has been inferred from a combination of astronomical observations which includes distance measurements of type la 
supernovae (SNe la) indicating that the expansion of the Universe is speeding up not slowing down (Perlmutter et al. 1999; Riess 
et al. 1998), CMB anisotropy data suggesting D.t - 1 (de Bernardis et al. 2000; Spergel el al. 2003), and clustering estimates 
providing ^ 0.3 (Calberg et al. 1996; Dekel, Burstein & White 1997). While the combination of the last two results impUes 
the existence of a smooth component of energy that contributes with ^ 2/3 of the critical density, the SNe la results require this 
component to have a negative pressure, which leads to a repulsive gravity. 

Despite the good observational evidence for the existence of these two forms of energy, it has never been shown that in fact 
they constitute two separate substances. In this concern, some authors have proposed the so-called Unified Dark Matter/Energy 
scenarios (UDME) or quartessence, i.e., models in which these two dark components are seen as different manifestations of 
a single fluid (see, for instance, Matos & Urena-Lopez 2000; Davidson, Karasik & Lederer 2001; Watterich 2002; Kasuya 
2001; Padmanabhan & Choudhury 2002). Among these theoretical proposals, an interesting attempt of unification was originally 
suggested by Kamenshchik et al. (2001) and developed by Bilic et al. (2002) and Bento et al. (2002), namely, the Chaplygin gas 
(Cg), an exotic fluid whose equation of state is given by 

Peg = -AIpI^, (1) 
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with a = I and A a positive constant. In actual fact, the above equation for or 1 constitutes a generalization of the original 
Chaplygin gas equation of state proposed by Bento et al. (2002). The idea of a dark-matter-energy unification from an equation 
of state hke Eq. (1) comes from the fact that the Cg can naturally interpolate between nonrelativistic matter (p = 0) and negative- 
pressure ip = -const.) dark energy regimes (see Bento et al. 2002; Alcaniz, Jain & Dev 2002 for details). 

Very recently, there has been a wave of interest in exploring theoretical (Bordemann & Hoppe, 1993; Hoppe 1993; Jackiw 
2000; Gonzalez-Diaz 2003a; 2003b; Kremer 2003; Khalatnikov 2003; Balakin et al. 2003; Bilic et al. 2003) and observational 
consequences of the Chaplygin gas, not only as a possibility of unification for dark matter/energy but also as a new candidate 
for dark energy only. These models have been tested for a number of cosmological data sets, including SNe la data (Fabris, 
Goncalves & de Souza, 2002; Cohstete Jr. et al. 2003; Avehno et al. 2003; Makler, de Oliveira & Waga 2003), statistical 
properties of gravitationally lensed quasars (Dev, Alcaniz & Jain 2003; Silva & Bertolami 2003), CMB measurements (Bento, 
Bertolami & Sen 2003a; 2003b; 2003c; Carturan & Finelli 2002; Amendola et al. 2003), age and angular size - redshift tests 
(Alcaniz, Jain & Dev 2002; Alcaniz & Lima 2003), measurements of X-ray luminosity of galaxy clusters (Cunha, Lima & 
Alcaniz 2003), future lensing and SNe la experiments (Avelino et al. 2003; Silva & Bertolami 2003; Sahni et al. 2003), as well 
as by observations of large scale structure (Multamaki, Manera & Gaztanaga 2003; Bilic et al. 2003). The present situation is 
somewhat controversial, with some tests indicating good agreement between observational data and the theoretical predictions 
of the model and others ruUng out the model as an actual possibiUty of description for our Universe (Sandvik et al. 2002; Bean 
& Dore 2003) (see, however, Bega et al. 2003). 

The aim of the current paper is to check the vahdity of such models with radio-selected gravitational lens statistics and also 
with a combination of gravitational lensing and SNe la data. To do so, we adopt the most recent radio source gravitational lensing 
sample, namely, the Cosmic All Sky Survey (CLASS) statistical data which consists of 8958 radio sources out of which 13 
sources are multiply imaged (Browne et al. 2002; Chae et al. 2002). Here, however, we work only with those multiply imaged 
sources whose image-splittings are known (or Ukely) to be caused by single galaxies, which reduces the total number of lenses to 
9. For the cosmological background we consider a flat scenario in which the generalized Cg together with the observed baryonic 
content are responsible for the dynamics of the present-day Universe (UDME or quartessence models). In our computations we 
adopt Sib = 0.04, in accordance with the latest measurements of the Hubble parameter (Freedman et al. 2002) and of the baryon 
density at nucleosynthesis (Buries, Nollett & Turner 2001). 

This paper is organized as follows. In Sect. 2 we present the distance formulae necessary to our analysis. In Sect. 3 we discuss 
the CLASS statistical sample, especially the observational criteria used as well as the restrictions adopted in our analysis. In Sect. 
4, we derive the corresponding limits on Cg scenarios from CLASS lensing statistics. We also examine the constraints obtained 
from the statistical combination of lensing data with recent SNe la observations and compare our constraints with others derived 
from independent analyses. Finally, in Sect. 5, we finish the paper by summarizing its main results. 



2. Basic equations 

By inserting Eq. (1) into the energy conservation law (u^T^J = 0) one finds. 



Peg = 



A + B 



or, equivalently. 



Peg = Peg, 



As + a-A 



,3(l+a) 



(2) 



(3) 



where pcg stands for the Cg energy density, the subscript o denotes present day quantities, Rit) is the cosmological scale factor, 
B = p^^" - A is a constant and A^ - Ajp^,^" is a quantity related to the sound speed of the Chaplygin gas today (vj = aAj). 
A fundamental quantity related to the observables here considered is the distance -redshift relation, given by 



X 



Roiio J(l+z)-i 



dx 



(4) 



y?-E(D.b, As, a, x) ' 

where x = ^ = (1 + z)~^ is a convenient integration variable, Clb stands for the baryonic matter density parameter, and the 



dimensionless function £(£1^, A„ a, x) is written as 



E = J%+(l-n,) 



A, -I- 



(1-A.) 



v.3(ff+l) 



1/2 



(5) 



For the lensing statistics developed in the next section, two concepts are of fundamental importance, namely, the angular 



diameter distance, Dis (zl, Zs ) 



R„ri (zL,zs ) 
a+zs) 



, between two objects, for example a lens at zl and a source (galaxy) at zs , 
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and the age-redshift relation, 

^ Ho Jo xE{Q.h,A„a,x)' 

^From the above equations, we note that UDME models reduce to the ACDM case for a - Q whereas the standard Einstein-de 
Sitter behavior is fully recovered for As -0 (see also Fabris et al. 2003; Avelino et al. 2003). 



3. Gravitational lensing statistics of the CLASS Sample 

Gravitational lensing directly probes the mass distribution in the Universe so that an investigation of lensing events of sources 
at high redshifts can provide important information about the global cosmological parameters and the structure of the Universe. 
The use of gravitational lensing statistics as a cosmological tool was first considered in detail by Turner, Ostriker & Gott (1984). 
Subsequently, it was realized that a comparison of theoretical lensing probabilities with gravitational lensing observations pro- 
vided an efficient constraint on the cosmological constant (Fukugita, Futamase & Kasai 1990; Turner 1990; Fukugita et al. 1992; 
Kochanek 1996) or more generally on the density and the equation of state of the dark energy component (Zhu 1998; 2000a; 
2000b; Sarbu, Rusin & Ma 2001; Chae et al. 2002; Huterer & Ma, 2003). However, the absence of an unbiased statistical sample 
of sources that is complete to within well-defined observational selection criteria and the uncertainties in the luminosity function 
(LF) of galaxies have seriously complicated the application of such methods. 

3.1. The CLASS Statistical Sample 

Recently, the CLASS collaboration ' reported the so far largest lensing sample suitable for statistical analysis, in which 13 out 
of the 8958 radio sources are multiply imaged (Myers et al. 2002; Browne et al. 2002). This sample is well defined through the 
following observational selection criteria (Myers et al. 2002; Chae 2002; Browne et al. 2002): (i) the spectral index between 
1.4 GHz and 5 GHz is flatter than -0.5, i.e., a > -0.5 with Sy oc v", where Sy is the flux density measured in milli-jansky; (ii) 
the total flux density of each source is > 20 mJy at 8.4 GHz; (iii) the total flux density of each source is > 5" = 30 mJy at 5 GHz; 
(iv) the image components in lens systems must have separations > 0.3 arcsec. The sources probed by CLASS at v = 5 GHz are 
well represented by a power-law differential number-flux density relation: \dN/dS \ oc (5/5°)'' with rj = 2.07 + 0.02 (1.97 + 0.14) 
for S > {< S^). The redshift distribution of unlensed sources in the sample is adequately described by a Gaussian model with 
a mean redshift z - 1.27 and dispersion of 0.95 (Chae 2002). Guided by the above information about (i) the number-flux density 
relation and (ii) the redshift distribution of unlensed sources, we simulate the unlensed radio sources (8945 in number) of the 
CLASS statistical sample using the Monte-Carlo technique (rejection method). 

In this paper, following Dev, Jain & Mahajan (2003), we work only with those multiply imaged sources whose image- 
splittings are known (or likely) to be caused by single galaxies. This means that our database is constituted by 9 lenses out of a 
sample of 8954 radio sources. 



3.2. Lensing Statistics 

We start our analysis by assuming the singular isothermal sphere (SIS) model for the lens mass distribution. As has been discussed 
elsewhere this assumption is a good approximation to the real mass distribution in galaxies (see, e.g.. Turner, Ostriker & Gott 
1984). In this case, the cross-section for lensing events is given by 

crsis = Ibn (-) (— ) , (8) 

c Dos 

where v represents the velocity dispersion and Dql, Dqs and Dis are, respectively, the angular diameter distances from the 
observer to the lens, from the observer to the source and between the lens and the source. By ignoring evolution of the number 
density of galaxies and assuming that the comoving number density is conserved, the differential probability of a lensing event 
can be expressed as 

cdt 

dr = «o(l -H Zifcrsis—dzL, (9) 

dZL 

where the quantity dt/dzi can be easily obtained from Eq. (7) and the present-day comoving number density of galaxies is 

n„ = I (f>(L)dL. (10) 
Jo 

' The Cosmic Lens All Sky Survey (CLASS): _http://wwvy.aoc.nrao.edu~smyers/class.htmI| 
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The differential optical depth of lensing in traversing dzt with angular separation between 
& Kasai 1990; Turner 1990; Fukugita et al. 1992): 

\2 



and (p + d(p is (Fukugita, Futamase 



-MdzL ^ F* {I +zl) 



dzidcp 
where the function F 



DmP, 



RoDos 
is defined as 



1 cdt 
R'odz[ 



X 



r/2 



r(a + i + p 



Dos 
Dls' 



X exp 



Dr 



D 



'LS 



d<p 



dZL, 



c 



4 , 
ff + - + 1 

r 



In Eq. ((T^, is the Schechter LF (Schechter 1976) given by 

L \" dL 
exp(-L/L,) — . 



<^(L)£/L = 0.| — 



(11) 



(12) 



(13) 



In order to relate the characteristic luminosity L, to the characteristic velocity dispersion v», we use the Faber- Jackson relation 
(Faber & Jackson 1976) for E/SO galaxies (L* oc vj), with j - A. For the analysis presented here we neglect the contribution of 
spirals as lenses because their velocity dispersion is small when compared to ellipticals. 

The two large-scale galaxy surveys, namely, the 2dFGRS ^ and the SDSS ^ have produced converging results on the total LF. 
The surveys determined the Schechter parameters for galaxies (all types) at z < 0.2. Chae (Chae 2002) has worked extensively on 
the information provided by these recent galaxy surveys to extract the local type-specific LFs. For our analysis here, we adopt the 
normalization-corrected Schechter parameters of the 2dFGRS survey (Folkes et al. 1999): a - -0.74, (p* - 0.82 x 10"^/z^Mpc~^, 
V* = 185 km/s and F* = 0.014. 

The normalized image angular separation distribution for a source at zs is obtained by integrating ^5"^ over zl- 



<yP_ 

d<p 



1 



r(zs) 



dh 



-dZL- 



(14) 



dzhdip 

The corrected (for magnification and selection effects) image separation distribution function for a single source at redshift 
Zs is given by (Kochanek 1996; Chiba & Yoshii 1999) 



P'(A0) = S 



y 

2A6 



r 

Jo 



Dos 
Dls' 



cdt 

x exp 

dZL 



Dos 
Dls' 



(1+zl)^ 



r(«-Hi + i) 

Similarly, the corrected lensing probability for a given source at redshift z is given by 

P'^TiZs) 



DqlDl 

RoDns 



1 



dZL- 



rdp 
J # 



Sd^. 



(15) 



(16) 



Here <p and AO are related as 



Afl 



and S is the magnification bias. This is taken into account because, as widely known. 



gravitational lensing causes a magnification of images and this transfers the lensed sources to higher flux density bins. In other 
words, the lensed sources are over-represented in a flux-limited sample. The magnification bias S(zs,Sv) increases the lensing 
probability significantly in a bin of total flux density {Sy) by a factor 



S(zs,S,) 



dN,,i>Sy) 



dSy 



r 



dN,,{>Sy/fl) 



dSy 



P(P) 



— dji. 
1^ 



(17) 



Here N7s{> Sy) is the intrinsic flux density relation for the source population at redshift zs- N,^{> Sy) gives the number of 
sources at redshift zs having flux greater than Sy. For the SIS model, the magnification probability distribution is p(fi) = ^/fr'. 
The minimum and maximum total magnifications and yu,,,,, ^ in Eq. ( I17> depend on the observational characteristics as well 
as on the lens model. For the SIS model, the minimum total magnification is ^ 2 and the maximum total magnification is 
yu„,„j = oo. The magnification bias S depends on the differential number-flux density relation \dNrj(> Sy)/dSy\. The differential 
number-flux relation needs to be known as a function of the source redshift. At present redshifts of only a few CLASS sources 
are known. We, therefore, ignore redshift dependence of the differential number-flux density relation. Following Chae (2002), 
we further ignore the dependence of the differential number-flux density relation on the spectral index of the source. 

An important selection criterion for the CLASS statistical sample is that the ratio of the flux densities of the fainter to the 
brighter images is > 0.1. Given such an observational limit,the minimum total magnification for double imaging for the 
adopted model of the lens is (Chae 2002): 

1 ^min 



= 2- 



(18) 



1 ^min 

Another selection criterion is that the image components in lens systems must have separations > 0.3 arcsec. We incorporate this 
selection criterion by setting the lower limit of AO in Eq. M6\ as 0.3 arcsec. 



2 The 2dF Galaxy RedshiftSurvey (2dfGRS): http ://msowww.anu.edu.au/2dFGRS/l 
^ Sloan Digital Sky S urvey : _http ://w w w. sdss .org/. 
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Fig. 1. Contour for 9 lensed radio sources in the parametric space A, - a. While the entire range of a is allowed, the parameter 
As is restricted to the interval 0.48 < < 0.66. 

4. Testing Cg Scenarios against Observations 

The expected number of lensed radio sources is Aliens = E P'i^ where P| is the lensing probability of the source and the sum is 
over the entire adopted sample. The expected number of lensed sources is thus a function of the parameters A , and a. We have 
done a grid search for those combinations (A,, a) by fixing Mens = 9. In Fig. 1 we show the contour for 9 lensed radio sources in 
the parametric space A , - a. As can be seen, while the entire range of a is allowed, the parameter A^ is tightly restricted to the 
interval 0.48 < A , < 0.66. This particular range for A , is not compatible with the one obtained from a SNe analysis involving 92 
events of the Supernova Cosmology Project and High-z supernova Search Team, i.e.. A , > 0.69 at 95% confidence level (Avelino 
et al. 2003) and is only marginally compatible with the limits from age + SNe performed by Makler, de Oliveira & Waga (2003). 
A comparison between the above interval with the one restricted by age estimates of high-z objects shows that the Mens test for 
the CLASS sample is compatible with the existence of the radio galaxy LBDS 53W091 (3.5 Gyr at z = 1.55) which implies 
A J > 0.52 but that it is not in accordance with the existence of the 4.0-Gyr-old radio galaxy 53W069 (at z = 1.43) and the 
2.0-Gyr-old quasar APM 08279+5255 (at z = 3.91) which requires, respectively. A, > 0.72 and A, > 0.82 (Alcaniz, Dev & Jain 
2003). The above interval from the Mens test is also not in agreement with the tight limits obtained by Silva & Bertolami (2003) 
from future SNe and lensing data, i.e., 0.75 < A., < 0.79 at la (or < 0.2). 
The Ukehhood function for lensing can be written as 



Here Ni is the observed number of multiple-imaged lensed radio sources and Nu is the number of unlensed sources in the 
adopted sample. The results of our analysis for the Cg model are displayed in Fig. 2a. The contours correspond to the 68.3% 
and 95.4% confidence level (cl) in the {As,a) plane. Although the entire range of a is permissible within the 68.3% confidence 
level, the parameter A, is constrained to be < 0.8 at 68.3% cl and < 0.9 at 95.4% cl. For this analysis the best fit model occurs 
for As - 0.36 and a - 0.35, which corresponds to an accelerating scenario with a deceleration parameter qo - -0.02 and a 
total expanding age of 7.33/i"' Gyr. Although not very restrictive, the constraints on the parameter Aj from the CLASS lensing 
sample of radio sources are more stringent than those obtained from the optical gravitational lensing surveys for quasars (Dev et 
al. 2003). However, it still limits the parameter a very weakly. 

4. 1. Joint analysis witii supernova data 

In order to perform a joint analysis with CLASS and SNe data sets, we first follow the conventional magnitude-redshift test 

(see, for example, Goliath et al. 2001; Dicus & Repko 2003; Padmanabhan & Choudhury 2003) and use the SNe la data set that 
corresponds to the primary fit C of Perlmutter et al. (1999) together with the highest redshift supernova observed so far, i.e, the 
I997fF at z = 1 .755 and effective magnitude m^^^ = 26.02 ± 0.34 (Benitez et al. 2002) and two newly discovered SNe la, namely, 




(19) 



1=1 



k=l 
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Fig. 2. a) Confidence regions in the plane A , - a arising from CLASS lensing statistics. Dashed and solid lines indicate contours 
of constant likelihood at 68.3% and 95.4% confidence levels, respectively, b) The likelihood contours in the a -Ag plane for the 
joint lensing + SNe la analysis described in the text. 

SN 2002dc at z = 0.475 and m'^f = 22.73 + 0.23 and SN 2002dd at z = 0.95 and m'^f = 24.68 + 0.2 (Blakeslee et al. 2003). We 
thus work with a total of 57 supernovae. The apparent magnitude of a given SNe is related to the luminosity distance di by the 
known relation m = M + 5 logDi, where = Hq cIl and M is given by the intercept obtained by fitting the low-redshift data 
set to m(z} = + 5 log(cz) (Hamuy et al. 1996). The value obtained is Al = -3.325 and confirms the results of Perlmutter et al. 
(1999). For the sake of completeness, we perform such SNe analysis. The best fit model occurs for As - 0.52 and a - -0.2 with 
a minimum value of = 64.31 (which corresponds to J = 1.17). When this magnitude-redshifttest is combined with CLASS 
lensing statistics, tighter constraints on the A , parameter can be obtained. As was shown, the index a is highly insensitive to SNe 
la data (Makler, de Oliveira & Waga 2003). Figure 2b shows the result of this joint analysis. For the combined analysis we 
\xsed.x%,gi = xlffe ~ ^l"^'' where I - £.iensl£l™a\ '^^e normalized likelihood for lenses. As can be seen, the limits on A., are more 
restrictive now than those imposed by the gravitational lensing statistics of the CLASS sample (Fig. 2a). Within 68.3% cl, the 
constraints on the parameter are as follows: 0.39 < Aj < 0.71 at 68.3% cl and 0.35 < A^ < 0.74 at 95.4% cl. In particular, the 
best fit model occurs for A^ = 0.58 and a = 0.5, corresponding to a 7.93/j~^-Gyr-old, accelerating Universe with a deceleration 
parameter = -0.33. 

5. conclusion 

A considerable amount of observational evidence suggests that the current evolution of our Universe is fully dominated by 
two dark components, the so-called dark matter and dark energy. The nature of these components, however, is a tantalizing 
mystery at present, and it is not even known if they constitute two separate substances. In this paper we have investigate some 
observational predictions of cosmologies driven by an exotic component named the generalized Chaplygin gas. These models 
constitute an interesting possibility of unification for dark matter/energy, where these two dark components are seen as different 
manifestations of a single fluid (UDME). We have investigated observational constraints from lensing statistics on spatially flat 
UDME scenarios. Since gravitational lensing statistics constitutes an independent way of constraining cosmological parameters 
we have used the most recent lensing data, namely, the Cosmic All Sky Survey (CLASS) sample to obtain the 68.3% and 95.4% 
confidence intervals on the parameters of the Cg equation of state. Our statistical analysis shows that the best fit scenario for these 
data occurs at As = 0.36 and a = 0.35. At 68.3% cl, parameter As is restricted to < 0.8 while the entire range of a is allowed. 
By considering the observed number of lensed radio galaxies we tightly constrain A, to the interval 0.48 < Aj < 0.66. From a 
joint analysis with SNe la data we obtain 0.35 < As < 0.74 at 95.4% cl with the best fit model occurring for As = 0.58 and 
a = 0.5, which corresponds to an accelerating scenario with a deceleration parameter qo = -0.33 and a total expanding age of 
7.93/i~' Gyr. As has been commented earlier, such results are only marginally consistent with those obtained from independent 
cosmological tests. It means that only with a more general analysis, possibly a joint investigation involving difi'erent classes of 
cosmological data, it will be possible to delimit the Aj - a plane more precisely, as weU as to test more properly the consistency 
of these scenarios as a viable possibility of unification for the dark matter and dark energy scenarios. 
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